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ABSTRACT 



We have analyzed the evolution of the spectral hardness parameter, E p k (the maximum of the 
vF v spectrum) as a function of fluence in gamma-ray bursts. We fit 41 pulses within 26 bursts 
with the trend reported by Liang & Kargatis (1996) which found that E p k decays exponentially 
with respect to photon fluence <fr(t). We also fit these pulses with a slight modification of this 
trend, where E p k decays linearly with energy fluence. In both cases, we found the set of 41 pulses 
to be consistent with the trend. For the latter trend, which we believe to be more physical, 
the distribution of the decay constant <I>o is roughly log-normal, where the mean of logio$o is 
1.75 ± 0.07 and the FWHM of logio'&o is 1.0 ± 0.1. Regarding an earlier reported invariance in 
$o among different pulses in a single burst, we found probabilities of 0.49 to 0.84 (depending on 
the test used) that such invariance would occur by coincidence, most likely due to the narrow 
distribution of $0 values among pulses. 

Subject headings: gamma-ray: bursts, observations 



Introduction 



The discovery of a gamma-ray burst optical counterpart with a measurable redshift seems 



to have shown that the sources are cosmological in origin ( Djorgovski et al. 1997 , Metzger 
et al. 1997). The observed fading multi- wavelength afterglows are so far consistent with the 



simple relativistic blastwave model (Meszaros & Rees 1997) which radiates via a synchrotron 
shock. However, the emission mechanism resulting in the prompt gamma rays remains a mystery. 
Studies of gamma-ray burst (GRB) spectral evolution have uncovered many trends which may 
be used to test possible emission mechanisms. In general, studies of GRB spectral evolution have 
focused on the "hardness" of bursts, measured either by the ratio between two detector channels 



or with more physical variables such as the spectral break or peak power energy E p k (Ford et 



al. 1995) which is the maximum of vF v , where v is photon energy and F v is the specific energy 



flux. Such hardness parameters were found to either follow a "hard-to-soft" trend (Norris et al 



1986| ), decreasing monotonically while the flux rises and falls, or to "track" the flux during GRB 



pulses ( Golcnctskii et al. 1983 , Kargatis et al. 1994 ). 

The discovery that E p k often decays exponentially in bright, long, smooth BATSE GRB 
pulses as a function of photon fluence $ (= f t ,Z FN(t')dt') (|Liang fc Kargatis 1996 , hereafter 
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LK96) provided a new constraint for emission models (Liang et al. 1997, Liang 1997, Daigne & 



Mochkovitch 1998 ). In their analysis, LK96 fit the function 



Epk(t) = E pk ( )e' 



-<£>(t)/<^ K 



(1) 



to 37 GRB pulses in 34 bursts. To interpret this empirical trend, they differentiated Eq. 1 to 
find 

- dE pk /dt = E pk Fn/^k ~ F E /^K (2 ) 



where F E = /^Xv^ E N ( E ) dE is the BATSE energy flux (see Eq. 1 of LK96). In this paper, 
we wished to avoid the assumption that E pk Fn ~ F E . To do this, we directly tested the trend 
— d(E p k)/dt = Fe/$o by integrating it to give us the function 



Ep k (t) = E pk(0) - £ (t)/$ 



(3) 



where £(t) (= J t ,~ , F E (t')dt') is the BATSE energy fluence. We emphasize that this is not a 
fundamentally different trend from the form used in LK96. 

The decay constant $>q K appeared to be invariant among pulses during some bursts analyzed 
in LK96, suggesting that individual pulses in a burst may originate in the same plasma. These 
discoveries coupled with the observed evolution of the spectral shape ( prider, Liang fc Preece 



1998a) suggested that saturated inverse Comptonization may be a viable mechanism during the 



gamma-ray active phase of bursts ( Liang et al. 1997 ), regardless of the distance scale (Liang 
1997|). 



Procedures 



To determine the evolution of GRB spectral shapes, we examined High Energy Resolution 
data collected from the BATSE Large- Area Detectors (LADs) and Spectroscopy Detectors (SDs) 
on board the Compton Gamma-Ray Observatory (Fishman et al. 1989). We began with the 



126 bursts which appear in Preece et al. 1998. These bursts were chosen for having a BATSE 
fluence (28-1800 keV) > 4 x 10~ 5 erg cm' 2 or a peak flux (50-300 keV on a 256-ms time scale) 
> 10 photon s _1 cm -2 . The counts from the detector most nearly normal to the line of sight of 
each burst (burst angle closest to 0) were background-subtracted and binned into time intervals 
each with a SNR of ~ 45 within the 28 keV to 1800 keV range. Such a SNR has been found to 



be necessary in time- resolved spectroscopy of BATSE gamma-ray bursts ( |Preece et al. 1998| ) . 

We deconvolved the gamma-ray spectra of each time interval using the Band et al. (1993) 
GRB function 



N E (E) = 



A ( E 

V 100 kcV 

A ' 



(a-|3)E 
100 kcV 



exp 

a-P 



(-£). («• 

exp(/3 - a) (joo^v)' 3 , (a • 



/3)E > E, 
/3)E <E, (4) 



where A is the amplitude (in photon s _1 cm~ 2 kev -1 ) and E = E pk /(2 + a). While LK96 
assumed that a and j3 were constant during the course of each burst, this has since been shown 
to be untrue with a larger data set. ( Crider et al. 199~7| ) . We thus left a and /3 as free parameters 
in our fits. 

At this point, we needed to select pulses within our bursts that we could use to test Equations 
1 and 3. Ideally, our pulses would not overlap other pulses and our method for choosing the time 
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bins associated with a pulse would not be biased. Unfortunately, by forcing our time bins to 
have a SNR ~ 45 so that spectra may be fit to them, much time resolution is lost. Pulses which 
would be easily separable at a higher time resolution become blurred together. In Figure ^, we 
show an example of what would likely be identified as two pulses in our coarse (SNR ~ 45) data. 
Below it, 64-ms count rate data for this same burst, obtained from the Compton Observatory 
Science Support Center (COSSC), is plotted. With higher time resolution, we see this bursts is 
composed of at least 4 distinct pulses. 

To avoid contaminating our sample with overlapping pulses (and to avoid biases introduced 
by a human in pulse selection), we used the COSSC 64-ms count rates and background fits to 
determine where each of our pulses began and ended. To do this, we developed an interactive 
IDL routine to fit the Norris et al. (1996) pulse profile to the individual pulses within these 
bursts. The pulse profile function for the count rate C(t) can be written 



where t max is the time of maximum count rate, tr r and o~d are the count rate rise and decay time 
constants, and v is the pulse "peakedness" . For an exponential rise and decay, v = 1. When 
v = 2 and er r = <7d this shape describes a Gaussian. 

With 64-ms resolution, we found that in many bursts pulses overlapped in a fashion making 
them too complex for us to fit individual pulses. Other bursts contained pulses which could 
be resolved, but none of their pulses lasted long enough to span at least 4 time bins with SNR 
~ 45. For 13 of the bursts, no processed 64-ms data was available. We discarded bursts which 
fell into any of these three categories. This left us with 26 bursts. This is comparable to the 34 
multi-pulse bursts used in LK96, which included several pulses which appear to be overlapping 
when confronted with the 64-ms data. To avoid overlap in our own analysis, we used only bins 
which were dominated by a single pulse (at least 70% of the counts from one pulse) . Within our 
26 bursts, we identified 41 regions composed of at least 4 time bins dominated by a single pulse. 
The time bins selected for each pulse were consecutive in all but two cases (BATSE triggers 451 
and 3290). The 64-ms data for each of these two bursts suggests that a short pulse occured near 
the middle of a longer pulse, which forced us to fit two separate regions with a single decay law. 

Our next step was to test the E p k-fluence relations (Equations 1 and 3) with each of the 
selected pulses. Our motivation for emphasizing the E p k-energy fluence relation (Eq. 3) as 
opposed the Ep^-photon fluence relation (Eq. 1) is that we believe that the former represents 
a more physical quantity. It is possible (perhaps even likely) that the observed BATSE photon 
fluence is a poor representation of the bolometric photon fluence. The BATSE LAD energy 
window was designed to contain the peak of GRB energy spectra, not the peak of the photon 
spectra. By using energy fluence in place of photon fluence, we can avoid the shakier assumption 
that the BATSE LAD photon flux is proportional to the bolometric photon flux. LK96 had 
attempted this but found that statistical errors in £ were too large to be useful. This was a 
result of their fixing a and /3 in the spectral fitting. When we fit the time-resolved spectra with 
variable a and f3, we obtained much smaller errors for £ , which made testing the E p k — £ relation 
possible. Nevertheless, we also fit Eq. 1 to our pulses in this paper both for historical reasons 
and as a test of our interpretation. 




(5) 
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3. Results 

We fit both the E p k — $ and the E p k — £ relations to our 41 clean pulses using FITEXY 



(Press et al. 1992). Table 1 summarizes the results for each of the pulses in our sample. The 
first column is the BATSE trigger number. The second column is the burst name, which is also 
the date the burst triggered in the format YYMMDD. The third column is the number of the 
LAD which was used for processing. The fourth column lists the t max from the Norris function 
fit to the pulse. The fifth column is the energy fluence within the bins selected for fitting in 
units MeV cm~ 2 . The sixth column gives the number of bins selected for fitting. The seventh 
column is the fitted value $q K for each pulse defined in Eq. 1. The eighth and ninth columns 
are the fitted values of $o and E pk (o) for each pulse as defined in Eq. 3. Of course, <&o will only 
equal <&q K if E p kFN = F E . Since the latter is not strictly true, we find that <E> ~ $o K - For 
completeness, we also show the plots of E p k versus £ and their fits in Figure ||. 

From the x 2 and the number of fluence bins for each decay fit, we calculated the probability 
Q of randomly getting a higher \ 2 by chance. Thus, Q >~ 0.5 represents very good fits, while 
Q ~ represents poor fits. The Q values from fits of Eq. 3 to our pulses appear in the tenth 
column of Table 1. If E p k does indeed cool linearly with £ in all pulses selected for fitting, then 
when plotting the cumulative distribution of Q values, we would expect 10% of the pulses to have 
a Q less that 0.1, 20% of the pulses to have a Q less than 0.2, and so on. Figure |^ shows the 
cumulative distribution of Q values for our pulses with acceptable fits. An excess of pulses with 
very high Q values would suggest a biased pulse selection process. A Kolmogorov-Smirnov test 
(P = 0.18) applied to our distribution of Q values suggests that the set of 41 pulses is not too 
biased and roughly follows the distribution we would expect if all of them are consistent with a 
linear decay of E p k with respect to energy fluence. 

By fitting the E p k — fluence law to the full observable duration of each pulse and not just the 
flux decay phase, we could characterize our pulses as "hard-to-soft" . None of our pulses required 
a "tracking" classification, though many of the ambiguous pulses excluded from this study could 
be "hard-to-soft" or "tracking". Three of the pulses in our sample (BATSE triggers 2316, 3491, 
and 3870) contain pulses with negative values of $o. However, all three of these pulses are still 
consistent with a positive value of <&o. We remind the reader here that large absolute values of 
$0 (like those in these three pulses) correspond to pulses with very little change in E p k, where 

= —g^r- » 0. In such cases, small statistical errors in translate to very large statistical 
errors in <J> . Even if all pulses decay monotonically from hard-to-soft, we should expect to see a 
few pulses to have negative values of 4>o- Since all of our pulses are consistent with a monotonic 
decay in E p k, we adopt the hypothesis that all pulses behave this way for the remainder of this 
paper and drop these three pulses from our sample to simplify our calculations. 



3.1. Distribution of $ 

The distribution of fitted $o values appears in Figure ^. It is roughly log-normal where the 
mean of logi $ is 1.75 ± 0.07 and the FWHM of logi $ is 1.0 ± 0.1. This distribution likely 
suffers some selection effects. This becomes obvious when one realizes that $ n ~ — J^f . Wc 
see that the smallest absolute value of $0 is limited by the minimum energy fluence which allows 
one to fit spectra (about 1 MeV cm~ 2 from Table 1) and the energy window of BATSE (max 
AE p k ~ 1870 keV). There are no such limitations on the high side of this distribution, since 
AE p k I can be arbitrarily small and A£ is only limited by nature. 
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3.2. Testing the Invariance of $o Among Pulses within Bursts 

LK96 reported that the decay constant $q K sometimes remains fixed from pulse to pulse 
within some bursts. Such behavior would hint at a regenerative source rather than a single 
catastrophic event (such as Mcszaros fc Rccs 1993| ). However, the intrinsically narrow distribution 
of decay constants mentioned above and the relatively large confidence regions for each pulse's 
value of <&o suggest that many bursts would appear to have an invariant decay constant merely 
by chance. 

As done earlier with a larger, but less reliable, set of pulses ( Crider, Liang fe Preece 1998b ), 
we calculated three statistics for each multi-pulse burst to test the invariance of the E p k — flucncc 
decay constant. We compared two of each bursts' M pulses at a time using the statistic 

x? _ [*.(.)- W (6) 

J a 1 4- n 1 

a *o(i) ^ a *o(j) 

and then distilled the comparisons within each burst into a single statistic to represent that 
burst. These statistics are defined in Table 2. Each is tailored for different null hypotheses. The 
statistic G\ tests if at least two pulses in a burst are similar (and thus "invariant"), while Gi tests 
if all the pulses have a similar decay constant. G3 tests for cither a single good pairing or several 
moderately close pairings. We believe that this last statistic is the most reasonable for testing our 
results since it does not require that all pulses decay at the same rate (as G2 does) but also does 
not discard information about multiple pulses repeating (as G\ does). Finally, we calculated a 
table of probabilities P(G,M) for our goodness-of-fit statistics G based on a simple Monte Carlo 
simulation. We created synthetic bursts with pulse decay parameters randomly sampled from 
the observed distributions of $0 and a<s> /&o- To avoid any bias that intrinsic invariances would 
have on these distributions, we created them using only one pulse from each burst. The sample 
of bursts in this study has fewer bursts than previous works, and hence has fewer bursts with 
more than one pulse. The three versions of the G statistic defined above are equivalent when only 
two pulses appear in a burst. Thus for this sample, with only 3 of the 9 multipulse bursts having 
more than 2 pulses, these statistics are nearly equivalent. The high probability that the observed 
repetitions occurred by coincidence leads us to conclude that pulse decays are not invariant from 
pulse to pulse within bursts. Instead, we suggest that the distribution of $0 values seen in all 
bursts is narrow enough that an apparent invariance of <&o is inevitable in some burst. We came 
to the same conclusion when examining <I>q K ( Crider, Liang fe Preece 1998b|) . 



Discussion 



Out of the 26 bursts to which we could fit a time-evolving Band GRB function, all contain 
at least one pulse consistent with a linear decay of Epk with respect to energy flucncc. Of the 



41 pulses in these bursts, all are consistent (Q > 0.001; Press et al. 1992) with this decay 
pattern. This is also true when we fit the LK96 exponential decay of E p k with respect to 
photon fluence. Besides LK96, other quantitative spectral evolution trends have been reported for 
GRBs. The averaged temporal and spectral evolution for 32 bright GRBs has been calculated in 
Fcnimorc 1998 . The averaged photon flux evolution can be described as both rising and decaying 
linearly with time. The hardness, as measured by E p k with a and (3 held fixed, also appears 
to decay linearly with time during the averaged burst (E p k = E pk (o)(l — t/t )). This is clearly 



not representative of all bursts since the evolution in bursts of E p k is often complex (Ford et al 



1995| , Liang fc Kargatis 1996| ). These trends possibly reflect the physics dictating the envelope 
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of emission. The fact that LK96 found the E p k-fluence trend in many mingled pulses may result 
from the fact that the burst envelope also evolves in this manner. 

Since the hardness of this envelope appears to decay more slowly than the hardness during 
the pulses we observe, we might not expect to see this trend in our pulses. However, the degree 
of confidence of E p k in our fits, coupled with the fact that energy fluence is often linear in time, 
makes the observations of many bursts possibly consistent (Q > 0.001) with this decay law. 
Testing the distribution of Q values as we did in Fig [|, we find a probability P=0.001 that the 
pulses are realizations of linear Epk-time trend, compared to P=0.18 for the linear Epk-energy 
fluence trend. While the linear E p k-time relation does not seem to describe individual pulses as 
well as the E p k-fluence relation, the results are not conclusive. 

More pulses are clearly needed if one is to discriminate between any two time-dependent 
spectral functions. One could simply wait for bursts to occur or for a more sensitive instrument 
to be built. However, it may be possible to increase the number of fittable pulses using the existing 
BATSE database. Fitting a time-dependent spectral function directly to higher time resolution 
data (or time-tagged event data) greatly reduces the number of required fit parameters. Another 
approach may be to analytically integrate the time-dependent spectral function and fit that to 
integrated spectra, as has been done by Ryde k Svensson (1998). By increasing the number of 
pulses, it will become possible to make more definitive statements about the evolution of prompt 
GRB emission and how it relates to the GRB afterglow. 
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pleasure to mention very useful discussions with Ed Fenimore, Charles Dermer, Markus Bottcher, 
and Rodney Hurley. This work is supported by NASA grant NAG5-3824. 
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Fig. 1. — The time history of BATSE trigger 543 seen both in (a) the lower time resolution 
which allowed us to fit a spectrum to each bin and (b) 64-ms resolution. Count rate is marked 
in each plot as a histogram and E p k, as determined from the data in the upper plot, is marked 
on both plots for convenience. While it is clear in the lower plot that there are at least 4 pulses 
within this burst, only two peaks in flux are evident in the upper plot. The fits of the Norris 
function (Eq. 5) to these pulses are plotted here as dotted lines. 

Fig. 2. — The E p k (crosses) and photon flux (histogram) evolution of 26 gamma-ray bursts with 
respect to energy fluence. Bins used to calculate the decay rate are marked with solid circles. 
The best fit decay law also appears. One sigma confidence bars for E p k and energy fluence are 
shown for each bin. 

Fig. 3. — Cumulative distribution of Q values for 41 pulses. In this plot, N is the fraction of 
pulses that have a Q value less than that of a certain pulse. If all of these GRB pulses cooled 
linearly as a function of energy fluence, one would expect the cumulative Q distribution to match 
the unit distribution (N=Q). The Kolmogorov-Smirnov test gives a probability P = 0.18 that 
this distribution is drawn from the unit distribution. From this, we conclude that our E p k — £ 
function (Eq. 3) adequately describes the pulses in this subset. 

Fig. 4. — This is the distribution of <&o values based on 41 pulses selected from 26 bursts. In this 
plot, we show both a histogram of values and the best-fit Gaussian distribution. The distribution 
of <f>o values is roughly log-normal where the mean of logio$ is 1-75 ± 0.07 and the FWHM of 
logi $ is 1.0 ±0.1. 
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Table 2: The probabilities of getting such good values for each of our three goodness-of-fit statistics G. These 
probabilities are high enough to suggest that any invariance seen in the data is purely coincidental. 

Definition of G a P(Gq( random ) < Ga( obscrvcd )) 
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